Density matrix of a single photon produced in parametric down conversion derived 
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We introduce an effective numerical method of density matrix determination of fiber coupled single 
photon generated in process of spontaneous parametric down conversion in type I non-collinear 
configuration. The presented theory has been successfully applied in case of source exploited in 
Phys. Rev. Lett. 99, 123601 (2007) to demonstrate the experimental characterization of spectral 
state of single photon. 
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INTRODUCTION 



The source of single photons is a prerequisite for im- 
plementation of linear optical quantum information pro- 
cessing schemes [Q, ||, ^, ^, quantum teleportation ||^ 
and quantum cryptography protocols |§ . Altough deter- 
ministic single photon sources has been recently devel- 
oped 10, H, |§, |lo| , most of the fundamental experiments 
still utilize the process of Spontaneous Parametric Down 
Conversion (SPDC). In this phenomena pairs of daughter 
photons are produced in a pumped nonlinear crystal and 
typically coupled into single mode fibers (SMFs). SPDC 
provides unique ease of shaping spectral and spatial prop- 
erties of generated nonclassical light. Detection of one 
photon, conventionally called signal, heralds the presence 
of the other, called idler |Q. Its statistical properties are 
satisfactory for post-selection type experiments as long as 
the pump power is low enough to make multiple pair gen- 
eration events negligible fl^. However, the idler photon 
is typically prepared in a mixed state ||l^ with respect to 
the temporal degree of freedom. Only by careful spectral 
filtering or precise adjustment of collecting setup one can 
produce pure wave packets jlj, |l^ [l^, Furthermore, 
coupling into SMFs decreases significantly the total num- 
ber of useful heralded photons. In a typical experimental 
scenario coupling only one photon out of the pair is signif- 
icantly more probable than coupling both photons. This 
comes at an expense of loosing sub-poissonian statistics 
of the counts but is beneficial for a certain test exper- 
iments or characterising vital characteristics of the 
source. 

We report here an effective numerical method of de- 
termination of the state of single photon generated in 
process of type I SPDC in configuration in which idler 
photon is coupled into SMF and the signal is disregarded, 
adopted from the results of jl^. This scheme has been 
previously used to demonstrate the method of character- 
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isation of a single photon state in Ref . ||l8| . The theoret- 
ical method is based on an observation that optical fiber 
define a relatively narrow range of directions that need to 
be included in the calculations. This justifies the paraxial 
approximation, which makes a substantial portion of the 
problem tractable analytically. In consequence we are 
able to derive the expression for the spectral density ma- 
trix of idler photon pi{ijJi,ijj[) propagating in SMF. Finally 
we present a numerical calculations of pi{uji,uj[) for a 
typical experimental settings. Both obtained results — 
highly simplified analytical assessment of the single pho- 
ton characteristics and an efficient method of calculating 
its full density matrix pi{ijji,u!'j) — can provide a valu- 
able insight in a future experiments with single photons 
obtained from SPDC. 



II. THEORETICAL MODEL 

Below we will develop a set of effective approximations 
for calculating the reduced density matrix pi{uJi,Lo'j) in a 
single arm of typical fiber coupled single-pass parametric 
down-conversion source. We start by writing down the 
full biphoton wavefunction in the perturbative approxi- 
mation. We describe the coupling of one of the photons 
into a SMF and we trace over the other photon. This 
way we obtain an involved, multidimensional integral for 
the reduced density matrix pi[ijji,u}[). At this stage we 
point out, that by judicious reordering of the integrals 
and performing a paraxial approximation this expression 
can be largely simplified and cast into a form which is 
easy to evaluate. 

Let us begin by describing the typical source depicted 
in Fig. |l|. It comprises the nonlinear crystal of length L, 
pumped by a beam of ultrashort pulses centered around 
frequency 2ujq. The elementary process taking place in 
nonlinear crystal is a decay of a pump photon into sig- 
nal and idler photons. We chose to parameterize the 
interacting waves by their frequencies Wp, Ws, Wj and 
transverse components of their wave vectors kp^, ksx, 
ki^. The energy conservation principle is expressed in 
this parametrization by an equation ujp = LOg -\- uji, while 
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Now we must take a certain spatial shape of the pump 
and the fiber spatial mode, preferably leading to ana- 
lytically integrable expressions. We assume there are no 
spatiotemporal correlation in the pump pulse and its am- 
plitude is a product of spectral Asp{uj) and spatial Up{'k±) 
parts: 



(3) 



We impose a spatial profile to be gaussian, corresponding 



to a beam of waist Wpi 



FIG. 1: The setup is comprised of x^^^ nonlinear crystal of 
length I/, pumped by a gaussian beam of Wp comprised of 
pules described by a spectral envelope function Aps. The 
idler photon is collected to a single mode fiber SMF, provided 
it is emitted into a spatial mode of characteristic width ui/at 
the crystal. 



the momentum conservation requires that the transverse 
components of wave vectors match kp^ = ks^ + k^j^. 
Note that for each of the waves once kj^ and uj are given 
kz is fixed. Therefore an ideal match of the longitudinal 
wave vectors is in general impossible. We choose a per- 
fect matching for pump photon of frequency ujp = 2loo 
propagating along the z axis (kp^ ~ 0) and signal 
and idler photons of equal frequencies ujs = uji = loq 
propagating symmetrically at the angle a with respect 
to z axis. The propagation directions of both gener- 
ated photons correspond to a transverse wave vectors 
IksOJ-l = IkioJ-l = wo/csina pointing in opposite direc- 
tions. Thus a phase matching, which we consider here, 
corresponds to fulfilling the following criterion: 

AA:^(kso_L,ti'o;kio_L,wo) = 

^p2(kso±+kioj_,2wo)-fcsz(kso±,wo)-fci^(kio_L,t^o) = 

(1) 

Let us write down an expresion for the probability am- 
plitude of generation two photons of frequencies Ws, i^i 
and transverse wave vectors ksj_ , k^^ jl^, |l^ : 

'I'(ksj_,tj5;k,;i,LJi) = 

fL/2 

N / dzAp(k,x+k,x,cc;,,+c^,;)e'^'^^^*'="'"='''-"'"-^', 

J-L/2 

(2) 

where N is normalization constant. The above formula 
expresses a fact that the probability amplitude of gener- 
ating pair of photons characterized by k^j. , oJs and k^^ , uji 
is proportional to the pump amplitude Ap{ksA_+'ki^,ujs + 
uJi). In turn the integration sums the contributions form 
the sHces of the crystal perpendicular to z axis, with ex- 
ponential term representing propagation phases. Indeed, 
Afcz (ks_L , Ws ; ki _L , ) z is the phase acquired by the pump 
while propagating to the sHce of interaction and by the 
photon pair while propagating towards the end of the 
crystal. 
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However, we still provide freedom of choice of the pump 
spectral amplitude Af^p{ijj). 

Next we find the wave function with respect to the 
fiber for an idler photon and with respect to the free 
space for the a signal photon. This can be done by pro- 
jecting i/'(ksj_, Ws; ki_L, Wj) on the fiber mode profile. We 
approximate the spatial mode coupled in SMF Mj(ksj_) 
by a gaussian function centered around the phase match- 
ing direction kio_L corresponding to a beam of width Wf 



Mj(k,;_L,LJi) 



— ^ exp I — ^ (k,^ - k,o±) 



(4) 



Without loss of generality we assume the collecting op- 
tics and the fibers are in the xz plane. Therefore the 
transverse wave vector of the idler photon is a function 
of its frequency oji in the direction of observation a is 
given by kioj. = —ijJiShi{a)x / c. Here a; is a unit vec- 
tor in the direction of x axis. Thus the biphoton wave 
function of idler photon inside SMF and signal propagat- 
ing in free space \ti(ksx, w^, uJi) is a partial overlap of the 
free space wave function ^'(ks±, cj^; 'ki^^LOi) and the fiber 
mode profile Ui{ksA_): 

rL/2 . 

*j(ks_L,u;s,Wi) =NAsp{uJs+>^i) / dz / d^k^^ 

J-L/2 J 

<(k,i,o.,)«p(k.,± +k,i)e'^'=^(''^-•'^=^'^-^"•)^ (5) 

In the next step let us write an expression for the den- 
sity matrix of the idler photon in the fiber pi{ijji,ijj[). It is 
obtained by taking a trace of the photon pair density ma- 
trix **(ks_L,Li;s;kij_,Wi)^'(k(,^,a;^;k^j^,w-) over the de- 
grees of freedom of the signal photon Ug and ks^ : 



= j dwgdksj. ^E'*(ks_L,a;s,Wi)\E'^(ks_L,a;s,^^■)■ 
(6) 

Thus a general structure of the spectral density matrix 
Pi{uji,uj[) after substitution of Eq. (||) into Eq. (^) is given 
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by: 



kioj^,k';^ — k'^Q^) to shorten the notation: 



Pi{uj.i, uj'i) (X J dzdz'dujsAsp{u;s + uJi)Asp{uis + w.')^ 

X J dks±dki±dk^j^u*UpU*u'^exp{iAk'^z' — iAkzz), 

(7) 

where the functions Up, Ui and Akz are taken for k^x, 
ki±, ojs, u!i and u^, m- and Afc^ refer to the respective 
quantities taken at k^x , k^_|_ , Ws , . A classical approach 
to evaluating this integral would be to first integrate an- 
alytically over the length of the crystal z and z' and then 
evaluate remaining 7 integrals numerically. However this 
would be extremely tedious computationally. Therefore 
we adopted another approach in order to alleviate the 
numerical effort. We first perform the paraxial approxi- 
mation of the phase mismatch exponent iAk'^z' — iAkzZ. 
Then the integrals in second line of Eq. (|^) can be evalu- 
ated analytically. Finally the remaining integrals over z, 
z' and ujs are computed numerically. 

The paraxial approximation of the phase mismatch ex- 
ponent iAk'^z' — iAkzZ is well justified as long as waists 
of the beam involved are much bigger than a wavelength. 
Mathematically this is equivalent to observing that in the 
integral Eq. (0) the range of k^x is in fact restricted by 
the aperture of the collection optics described by Uj;(ksx) 
to kix — kioj_ ± as seen in Eq. (^. Then the 

range of perpendicular component of the signal wave vec- 
tor ksx is also restricted because of finite wave vector 
spread of the pump. Direct calculation reveals the range 

of approximately |Q: k^x — k^ox ± y^l/w^ + l/w^, 
where k^ox = sin(Q;)i/c. Those two facts enable us to 
approximate the phase mismatch Akz{}^.s^_,^s]'ki^^^i) 
around the center of both ranges up to second order in 
ksx and k^x: 

AA:^(ks_L,Ws;kj_L,cjj) := 
kpz{ks± + k,x,t^s - fcs^(ksj_,Ws) - kiz{\<ii±,uji) « 

(8) 



Upu'pU*u[eyi\i{iAk'^z' — iAkzz) 



exp (-/«M2K + Mik; + Mo) (9) 

where Mq, Mi and M2 depend on frequencies Ug, uji, uj[ 
and z, z' . Their explicit form is given in Appendix |^. 
Now the integral in Eq. (|^) can be partially evaluated 
and the density matrix pi{uji,uj'^ is found to be: 



Pi{uji,uj'i) = \JVf / du!sA* Jujs +u!i)Asp{LUs 



L/2 L/2 



dz / dz' 



-L/2 -L/2 



, expMo / 1„,T, ,-1,., 
' exp — Mf Mo ^Mi 

Vdet M2 4 ^ 2 



(10) 



III. NUMERICAL RESULTS 



In order to proceed we fix the spectral shape of pump 
pulse Asp{uj). We assume it be gaussian of duration Tp-. 



Vt'p 



exp 



2 ^ 



(11) 



In typical experimental scenario photons are spectrally 
filtered in order to ameliorate its characteristics or reduce 
a noise. We can model spectral filtering effects multiply- 
ing the density matrix by the the amplitude transition 
functions A{u}'): 



Pi{uj,Lo') ~* A{uj)A{uj')pi{uj,uj') 



(12) 



We approximate amplitude transmission of spectral filter 
by a Gaussian function defined in the following way: 



A(c.)^cxp(-21n2 ^"-"f 



(13) 



(k. 



where for sake of brevity we introduced k 
ksO-Ljkix — kioj.)"^ to denote the four component vector 
of deviations form the directions of phase matching. The 
expansion coefficients AA;2(ks_L, k^x, o^A Di(a;s, 0;^) 
and 'D2i^^s,'-^i) are given in the Appendix |A[ Note that 
for each pair of frequencies uis and Ui we have in gen- 
eral different expansion coefficients. This means we keep 
dispersion relations exact. This is motivated by the fact, 
that the spectrum of the single photons can easily span 
lOOnm or more, which makes expansion of the k vector 
as a function of frequency inaccurate. 



With Taylor expansion given by Eq. (Al) at hand, we 
can perform inner integrals in the expression for density 
matrix Eq. (0). Let us rewrite the integrant of Eq. (0) 
using a vector of deviations k = 0<-s± — kso± , k^x — 



Here we have calculated the spectral density matrix 
of photon generated within the typical experimental sce- 
nario, where the nonlinear crystal is pumped by a gaus- 
sian pulses Eq. (11). Note that for real gaussian pulses 
density matrix is real. Fig. || shows the spectral den- 
sity matrices p{lu,uj') of photons generated in L = 1mm 
long BBO crystal cut for 9c ~ 30°, pumped using gaus- 
sian pulsed beam of Wp = 100/^m waist and duration 
Tp = lOOfs. The collecting optics and the single mode 
fiber were set to couple the gaussian mode of waist equal 
Wf ~ 200/im. The spectral filters of cr = 20nm were 
applied. Panel a) shows density matrix when the obser- 
vation direction coincides with the phase matching di- 
rection a = 2.2°, see Fig. 0. In panel b) the observation 
direction is set to a = 3.2°. 
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(a) 



(b) 
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FIG. 2: The spectral density matrix p{ijj,Lo') of photons trans- 
mitted through 20nm spectral filter. The plots were generated 
for L — 1mm long BBO crystal cut for 6c — 30°, pumped 
using gaussian beam of Wp = 100/im waist. The pulse dura- 
tion was set to rp = lOOfs. The collecting optics and the 
single mode fiber was set to couple the gaussian mode of 
Wf = 100/im. The observation angle was set to a) q = 2.2° 
(direction of perfect phase matching) and h) a = 3.2°. 



and 



where we denote the blocks as: 



and 



dkfj,r,dk 
d'^Ak, 



dk 



d^Ak, 



(A3) 



(A4) 



(A5) 



In the above indexes [i,v = s,i refer to signal and idler 
photons. 



IV. CONCLUSIONS 

We have derived a simple method of predicting the 
spectral density matrix piiiOi^uj'^ of a photon generated 
in one arm of pulsed SPDC. The model allows to re- 
duce dramatically the computational effort while retain- 
ing the accuracy of results. Those results, by the virtue 
of their simplicity, may become very useful when engi- 
neering photon sources to meed a certain needs for the 
spectral characteristics. 

The presented model has been successfully used to pre- 
dict the measurement outcome in Ref. Q . Note that 
the model Eq. (0) may be used to predict the spectral 
density matrix in case of any type of spectral pumping. 
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APPENDIX A: PARAXIAL APPROXIMATION 

The standard paraxial approximation is given by: 

AA:^(ks_L,Ws;ki_L,Wj) w Afc^(k30_L, Ws; kjo±, Wi)-!- 

+ Vt{(uJs,'^i) ■ K + K'^'D2{^s,UJi)n, (Al) 

where for sake of brevity we introduced k = (ks± — 
ksO-Ljkj^ ~ ^io±)'^ to denote the four component vec- 
tor of deviations form the directions of phase matching. 
Akz{'ksO±,i^s','^io±,^i) is spatially constant term while: 



Di(ws, Wi) 



(A2) 



APPENDIX B: DEFINITIONS OF DENSITY 
MATRIX COEFFICIENTS 

The formula for spectral density matrix can be written 
down in the very compact form: 



-L/2 



thanks to the following definitions: 

/2l2 
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©2 iz 


'd' 



Ml 
Mo 



-(Ao + Ao')Wp -izds- 
— AoWp — izdi 
-Ao'wl + iz'd'. 



(B2) 
(B3) 



w: 



^ (Ao^ + Ao'^) - izAkz{'ksOj_,uJs;'kto±,L^i) 



-iz'AA:^(kso±,t^s;k-| 



(B4) 



We use IL2 and O2 to denote two dimensional iden- 
tity and zero matrices, furthermore Aq = k^oj. + k.joj_, 
Aq' = kso_L + k-o± ^iid all primed Taylor series coeffi- 
cients components d^^ and d^ are taken for the primed 
idler frequency wj. 
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